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ABSTRACT 

o 
o 

*vj I We examine the reprocessing of X-ray radiation from compact objects by accretion 

disks when the X-ray emission from the star is highly beamed. The reprocessed flux 






D ' for various degrees of beaming and inclinations of the beam axis with respect to 

the disk is determined. We find that, in the case where the beam is produced by a 



> 



non-relativistic object, the intensity of the emitted spectrum is highly suppressed if 
the beam is pointing away from the disk. However, for beams produced by compact 
objects, general relativistic effects cause only a small reduction in the reradiated flux 
OO I even for very narrow beams oriented perpendicularly to the disk. This is especially 

—^ I relevant in constraining models for the anomalous X-ray pulsars, whose X-ray emission 

0^ ' is highly beamed. We further discuss other factors that can influence the emission 

^^ ' from disks around neutron stars. 
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C^ ■ 1. Introduction 



x 



Reprocessing of X rays by accretion disks is a common phenomenon in astrophysics. For 
Vh ' example, reprocessed radiation is thought to be responsible for the optical emission observed in 

X-ray binaries (e.g. van Paradijs & McClintock 1994; de Jong, van Paradijs & Augusteijn 1996) 
and the ratio between optical and X-ray flux has been connected with properties of the system 
itself, such as the orbital period (van Paradijs & McClintock 1994). 

Another class of astrophysical objects for which reprocessing of radiation by an accretion 
disk could be potentially important is that of the Anomalous X-ray Pulsars (AXPs). The 
luminosity of these objects, Lx ~ iO'^^~^^ erg s~^, is far larger than their rotational energy loss, 
\E\ = A-K^IP/P^ ~ 10^^'^ erg s^^, and therefore their X-ray emission cannot be rotation-powered. 
Competing models to explain the origin of their X-ray luminosity invoke either the existence of 
neutron stars (NSs) with very large magnetic fields, B ~ 10^^~^^ G (e.g. Duncan & Thompson 
1992; Thompson & Duncan 1996; Heyl & Hernquist 1998a, b) or accretion from either a companion 
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(Mereghetti & Stella 1995), or a disk; where the latter situation could have resulted from the 
debris of a disrupted high mass companion (van Paradijs, Taam &: van den Heuvel 1995; Ghosh, 
Angelini & White 1997) or from material falling back after the supernova explosion (Corbet et al. 
1995; Chatterjee, Hernquist & Narayan 2000; Alpar 1999, 2000; Marsden et al. 1999). So far, no 
evidence for a companion has been found for any AXP, and the debate on whether the AXPs are 
isolated neutron stars powered by internal energy or NSs accreting from a disk of some sort is still 
open. 

The strongest constraint on the presence of a disk is given by its optical and longer wavelength 
emission which, at the luminosities typical for AXPs, is dominated by reradiation. Phase 
observations of the X-ray luminosity show a high degree of pulsation, which cannot be accounted 
for by an isotropic X-ray emission, but require some degree of beaming (DeDeo, Psaltis & 
Narayan 2000; Perna, Heyl & Hernquist 2000). However, computations of the reprocessed flux 
for disks around AXPs (Perna, Hernquist & Narayan 2000; Hulleman et al. 2000) have so far 
always assumed isotropic emission from the star. In this Letter, we consider the problem of disks 
irradiated by beamed radiation. For radiation produced by a compact object, such as a neutron 
star, the general relativistic effect of light bending is taken into account in transforming the beam 
pattern formed within the atmosphere of the neutron star into that which is seen by an observer 
at infinity. We show that the inclusion of these relativistic corrections has a significant influence 
on the reradiated spectrum when the beam points away from the disk. We finally discuss, more 
generally, other effects that can influence the spectrum of the reradiated flux, particularly in the 
context of the types of disks expected in accretion models for AXPs. 



2. Model 

2.1. Beaming pattern 

The radiation pattern emerging from accretion into the atmosphere of a neutron star (NS) 
has been computed in detail by Nagel (1981) and Meszaros &: Nagel (1985). At low accretion 
rates, a good approximation is given by the function 

/(<5)=/ocos"(5), (1) 

with n ~ 2 — 3, where 6 is the angle that the emitted photons make with the normal to the local 
stellar surface. What we need to find is the pattern of radiation that an observer at infinity will 
detect. Due to light bending effects, a photon emitted at an angle 5 with respect to the normal to 
the surface comes from a colatitude 6 on the star given by (e.g. Page 1995) 




e{6) = xdu / d(l-^]{-^] - (1 - 2u)u^x^ , (2) 



where x = sin 6, R is the radius of the star, and Rs = 2GM/c^ its Schwarzchild radius. If the 
beam pattern formed in the atmosphere of the NS is I (6), then an observer at infinity will see a 



beam pattern given by 

Given a colatitude 0* on the star, the corresponding angle (5* at the emission point is computed 
by numerically solving the equation 9^, — 6(6) = 0. The intensity at that point is then found from 
Equ. (|l]), and converted to the observed radiation emitted at ^* with the transformation given in 
Equ. i). 

Figure 1 shows I{0) for a star of mass M = IAMq and different choices of radii. The more 
relativistic a star is (i.e. smaller R/Rg), the more spread out the radiation pattern becomes. 

Now, let a be the angle that the axis of the beam makes with the plane of the disk; the axis is 
then described by the unit vector O = (cosacosi?i>, cosasini?!), sina), where (p = fit, with fi being 
the rotation rate of the star. Then, let P = (sin/3, 0, cos /3) be a point along a direction making an 
angle /? with the perpendicular to the disk in the upper plane. While the star rotates, the angle 
formed between the beam axis and the direction along P is 

9up{(p, P;a) = arccos(sin/3cosacosi?i) + cos/3sina) . (4) 

Due to the symmetry of the problem, two beams are expected from the NS, centered around the 
two poles. If the axis of the upper beam forms an angle a with the plane of the disk, the axis of 
the lower beam is at an angle —a, and, while the star rotates, it makes an angle 

^down(';^, (3; a) = arccos(sin (3 cos a cos cp — cos /3 sin a) . (5) 

with the direction P defined above. 

Figure 2 shows the intensity profiles averaged over a rotation period of the star, as a function 
of the angle f3 of the observation point with respect to the normal to the disk, 

(/(/?; a)) = ^ r d<P {/[0up(0, /?; a)] + /[0down(0, /?; a)]} . (6) 



2.2. Temperature profile and emitted spectrum 

Following Vrtilek et al. (1990), the computation of the temperature profile is obtained by 
combining two equations for the disk structure. The equation for hydrostatic equilibrium in the 
z-direction is 

-^T = a;2/i2 , (7) 

where a; is the Keplerian rotation rate, u; = yjGM/r''^, h is the half-thickness of the disk, k the 
Boltzmann constant, jjl the mean molecular weight, and rup the proton mass. The second equation 
for the disk structure states that the radiation emitted by a surface element, cTjj, should be 



equal to the sum of the energy released locally and the absorbed radiation (Cunningham 1976; 
Pacharintanakul & Katz 1980) 



where 



^^--^^0^^^' («) 
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gives the effective temperature for a disk which is not irradiated (Shakura & Sunyaev 1973). In 
Equation ^, /i is the absorbed fraction of the radiation impinging on the surface (the quantity 
1 — /i is called the "albedo" of the disk), for which we assume /i = 0.5 as in Vrtilek at al. unless 
otherwise stated. The X-ray luminosity is given by 

L^{P;a) = Ll {I{p;a)) , (10) 

with (/(/3; a)) given by Equation (^), and normalized so that J^ (/(/?; a)) sin/? dp = 1. L^. is the 
luminosity corresponding to isotropic emission, and it can be related to M via 

L.^A^. (U) 

where a value of /2 = 0.5 corresponds to the assumption that the X-ray luminosity produced in 
the disk is radiated mainly in a direction perpendicular to the orbital plane, and that only the 
X-rays directly from the NS can impinge on the disk. Finally, the angle /? in Equation ( p!o| ) is 
related to the half-thickness of the disk at position r by 



p = arctan 

2 



(12) 



Note that the vertical structure of the disk is assumed to be isothermal. 



Equations (0) - (jT^) are solved numerically for the temperature profile T{r) , and the flux to 
the observer is obtained by integrating the local emissivity B^(T) (assumed black body) over the 
entire surface of the disk: 

1 /rmax /•2n 

F = -^ / drr / d(j) By{T) x [cos(/)sin7(r) sini -|- cos7(r) cosi] . (13) 

U Jrinin JO 

Here i is the angle between the line of sight and the normal to the disk at r = 0, and 
7(r) = dh{r)/dr is the tilt angle of the surface of the disk at position r. In our calculations we 
assume i = 60°. The inner edge of the disk, rmim is taken at the magnetospheric radius Rm: for 
a typical luminosity Lx ~ 10^^ ergs s~^ and a magnetic field i? ~ 8 x lO^^G, as required by the 
accretion model of Chatterjee et al., one has Rm ~ 2 x 10^ cm. The outer edge, .Rout, is determined 
according to the solution given by Cannizzo (1993); for a disk of about ~ 10^ years, as typical for 



the AXPs, it is0 -Rout ~ 5 x 10^^ cm. We show our results for an X-ray luminosity Lx = 10^^ ergs 
s-i. 

Figure 3 shows the emitted flux for the exponents n = 2 and n = 3 of the beaming function, 
and for the limiting inclination angles of the beam axis a = 0° and a = 90°. Each case is made 
for R/Rs — > oo (non-relativistic case), and for R/Rg = 3. The emission corresponding to the 
isotropic case is shown for comparison. The figure illustrates how, in the non-relativistic case, the 
reprocessed flux is suppressed by several orders of magnitude if the beam axis is pointing away 
from the disk. However, in the case of radiation produced by a highly relativistic object such as 
a neutron star, even a beam perpendicular to the disk (a = 90°) gives a significant contribution 
to the spectrum. This reduction in the difference between spectra produced by beams pointing 
parallel and perpendicularly to the disk is due to the same general relativistic effects that cause a 
reduction in the pulsation amplitude of the radiation from the neutron star itself (Page 1995). 



3. Discussion 

We have considered the effects that beaming of the X-ray luminosity impinging on a disk has 
on the flux that is reradiated by the disk itself. We have found that, due to general relativistic 
effects of light deflection, even narrow beams pointing away from the disk still give a significant 
contribution to the impinging flux and on the consequent intensity of the reprocessed radiation. 

The problem we have considered is particularly important for the AXPs, which, due to the 
high pulsation amplitude of their flux, must produce a highly beamed X-ray luminosity. The 
reradiated flux gives the dominant contribution at long wavelengths; this is produced in the bulk 
of the disk, whereas the optical emission is generated in its innermost part. If the inner edge 
of the disk were truncated at a radius larger than the magnetospheric radius (which could for 
example happen if the inner region of the disk made a transition to a different state), then the 
optical emission would be highly suppressed. This is shown in panel (a) of Figure 4, where the 
spectrum is computed for different choices of rmin. Constraints from optical observations have 
already been made in a few cases. HuUeman et al. have used I and R band observations to 
rule out accretion scenarios for the AXP 1E2259. However, if this object had an accretion disk 
truncated at rmin ~ WRm, the observational constraints would not be violated any longer^. 

On the other hand, we have shown that the longer wavelength emission is hardly suppressed, 
even if the X-ray emission is highly beamed and the beam axis points away from the disk. A 



"^We also assume that the disk has a high opacity, as it is hkely to be made of heavy elements if formed from the 
debris of a supernova explosion. 

^Note that these constraints are given assuming an inclination angle of 60° between the line of sight and the 
normal to the disk. Clearly, a higher inclination would not violate the observational constraints while allowing a 
larger rmin- 



reduction of the disk to 10% of its size (panel (b) in Figure 4), or an increase in the disk albedo 
from 50% to 95% (panel (c) of Figure 4) would not reduce the peak flux by more than an order 
of magnitude. Therefore, definite constraints on the nature of the AXPs can only be made from 
multiwavelength observations. 

We are indebted to Ramesh Narayan and Jeremy Heyl for stimulating discussions. 
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Fig. 1. — Beaming pattern as it would appear to an observer at infinity for a NS of mass M = IAMq 
and different values of its radius. The emitted pattern at the NS surface is I{6) = cos^ 6. Light 
deflection effects make this pattern much broader. 
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Fig. 2. — Average intensity over a period of the star as a function of the angle with the normal to 
the disk. Here R/Rg = 3, and in both panels the various curves correspond to different choices of 
the inclination angle a of the beam axis with the plane of the disk: a = 90° (solid line), a = 60° 
(dotted line), a = 30° (dashed line), a = 0° (dotted-dashed line). The intensity is normalized to 
the value of isotropic luminosity. 
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Fig. 3. — Emitted spectrum from a disk of age about 10^ yr, illuminated by an X-ray luminosity 
L^ = 10^^ ergs sec~^. In both panels, the various curves correspond to the following cases: isotropic 
luminosity (solid line), R/Rg = 3 and a = 90° (dotted line), R/Rg —>■ oo and a = 90° (dashed 
line), R/Rg = 3 and a = 0° (long-dashed line), R/Rg -^ oo and a = 0° (dotted-dashed line). 
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Fig. 4. — The effect on the emitted flux of reducing the inner [panel (a)] or outer [panel (b)] edge 
of the disk, or having a larger albedo [panel (c)] is shown for L^ = 10^^ ergs sec~^ and isotropic 
X-ray emission. 



